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Abstract 

The river-resident Salmo salar ("smablank") has been isolated from other 
Atlantic salmon populations for 9,500 years in upper River Namsen, Norway. 
This is the only European Atlantic salmon population accomplishing its entire 
life cycle in a river. Hydropower development during the last six decades has 
introduced movement barriers and changed more than 50% of the river habitat 
to lentic conditions. Based on microsatellites and SNPs, genetic variation within 
smablank was only about 50% of that in the anadromous Atlantic salmon 
within the same river. The genetic differentiation (_F ST ) between smablank and 
the anadromous population was 0.24. This is similar to the differentiation 
between anadromous Atlantic salmon in Europe and North America. Microsat- 
ellite analyses identified three genetic subpopulations within smablank, each 
with an effective population size N e of a few hundred individuals. There was no 
evidence of reduced heterozygosity and allelic richness in contemporary samples 
(2005-2008) compared with historical samples (1955-56 and 1978-79). How- 
ever, there was a reduction in genetic differentiation between sampling localities 
over time. SNP data supported the differentiation of smablank into subpopula- 
tions and revealed downstream asymmetric gene flow between subpopulations. 
In spite of this, genetic variation was not higher in the lower than in the upper 
areas. The meta-population structure of smablank probably maintains genetic 
variation better than one panmictic population would do, as long as gene flow 
among subpopulations is maintained. Smablank is a unique endemic island 
population of Atlantic salmon. It is in a precarious situation due to a variety of 
anthropogenic impacts on its restricted habitat area. Thus, maintaining popula- 
tion size and avoiding further habitat fragmentation are important. 



Introduction 

The understanding of genetic and structural processes in 
isolated populations is a central topic in conservation 
biology and genetics (e.g., Groom et al. 2006). During the 
history of life on Earth, geological processes have been 
important both in creating barriers between populations 
and in merging formerly isolated populations. Terrestrial 
and marine habitats are relatively continuous and may 
provide opportunities for gene flow among adjacent pop- 
ulations, in contrast to freshwater bodies, which are often 
discontinuous, and in many ways similar to islands. Many 
freshwater species therefore constitute a number of popu- 



lations, which are isolated within restricted geographical 
areas (e.g., Bronmark and Hansson 2005). This low 
genetic connectivity has resulted in a large number of 
freshwater species (e.g., about 50% of teleost fishes; 
Wootton 1998) relative to the small amount of freshwater 
on Earth (about 2% of available water). 

Fresh water bodies at higher latitudes harbor relatively 
few fish species. Since the last deglaciation, and the subse- 
quent isostatic rebound, topography and climate have 
restricted immigration opportunities, resulting in species- 
poor fish communities. One major group entering fresh- 
water systems in previously glaciated areas of both North 
America and Eurasia was the anadromous salmonids 
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(e.g., Bronmark and Hansson 2005). In northwestern 
Europe, three salmonid fish species increased their distri- 
bution area in early postglacial times around 11,000 years 
ago, colonizing watercourses as they became accessible. 
Arctic charr (Salvelinus alpinus [L.]) and brown trout 
(Salmo trutta L.) readily established freshwater-resident 
populations, while Atlantic salmon (Salmo salar L.) rarely 
did so. Atlantic salmon normally maintain an anadro- 
mous life cycle, as all females smoltify (i.e., undergo phys- 
iologic and morphological changes to facilitate life in 
seawater) and migrate from the nursery stream to a more 
productive feeding environment in the sea, before return- 
ing to their nursery stream to spawn as adults. Some 
males remain in the river and mature sexually at a small 
body size (e.g., Webb et al. 2007; Jonsson and Jonsson 
2011). 

A few Atlantic salmon populations are freshwater sta- 
tionary and use a freshwater lake for feeding and growing 
instead of migrating to the sea. This type of nonanadrom- 
ous life cycle is common in parts of North America 
(Webb et al. 2007; Jonsson and Jonsson 2011), while only 
few populations are known in Europe. Eight localities 
with nonanadromous populations utilizing river-lake sys- 
tems are known in Russia, one in Finland (Lake Saimaa 
area), one in Sweden (Lake Vanern), and two in Norway 
(River Otra/Lake Byglandsfjord and River Nidelva/Lake 
Nelaug, the latter being extinct) (Dahl 1928; Kazakov 
1992; Nilsson et al. 2001; Barlaup et al. 2005; Siiisa et al. 
2005; Jonsson and Jonsson 2011). A completely river-resi- 
dent life cycle is even more exceptional. Except water- 
sheds in Newfoundland, Canada (Gibson et al. 1996; 
Webb et al. 2007), it is only known from River Namsen, 
Norway (Berg 1953, 1985). River-resident populations 
may arise when both sexes mature at the presmolt stage. 
Presmolt maturation is common in Atlantic salmon 
males, but rare in females (Hindar and Nordland 1989). 
Hence, river-resident populations have likely arisen from 
a low number of females and been subject to a strong 
founder effect. Being river resident also implies that all 
age-groups share a restricted area in fresh water. Anadro- 
mous salmon have the potential for larger population 
sizes because they utilize the sea during part of the life 
cycle. Anadromous salmon may also receive immigrants 
from other rivers (Jonsson and Jonsson 2011). Conse- 
quently, river-resident populations are expected to lose 
genetic variation at a higher rate. 

The river-resident Atlantic salmon in River Namsen, 
Norway (Fig. 1; with the colloquial name "smablank," 
which is used hereafter) constitutes a unique entity 
among European Atlantic salmon. Other freshwater-resi- 
dent populations in Europe migrate to lakes, they have all 
been depleted, and most are presently maintained by 
stocking programmes (Ozerov et al. 2010; Jonsson and 




Figure 1. Live wild specimens of the river-resident Atlantic salmon 
(Salmo salar) "smablank" from River Namsen. Above is an adult 
female, veteran spawner, approximate body length and mass: 
180 mm, 70 g. Below is a smaller adult male. Photograph: Per H. 
Olsen, NTNU, Trondheim, Norway. 

Jonsson 2011). The smablank in River Namsen constitutes 
an island population, which has been isolated for approx- 
imately 9,500 years (cf. Frankham 1997). Earlier analyses 
have demonstrated that smablank exhibits a lower genetic 
variation than the anadromous Atlantic salmon in the 
lower part of the river (Stahl and Hindar 1988; Vuorinen 
and Berg 1989; Bourret et al. 2013). Allozymes analyses 
have also shown genetic differences between smablank 
populations (Vuorinen and Berg 1989). At the same time, 
there was no evidence of genetic signatures of offspring of 
anadromous Namsen salmon, which were released as fry 
in smablank territory during 1950-1976. By applying 
more fine-meshed methods (microsatellites and SNPs) to 
analyze samples from a number of localities within the 
restricted distribution area, and also to analyze samples 
collected over the last six decades, we aim to understand 
population viability and resilience of this unique Atlantic 
salmon population. 

The hydromorphological features of the river sections 
inhabited by smablank indicate an asymmetric gene flow 
from the upper to the lower regions (Hanfling and Weet- 
man 2006). Fish moving upstream will face more obsta- 
cles and barriers than fish moving downstream (Kawecki 
and Holt 2002). Asymmetric gene flow is a general phe- 
nomenon for many organisms living in habitats with pre- 
dominantly unidirectional water currents (Pollux et al. 
2009; Pringle et al. 2011), and indeed in most species 
with a source-sink metapopulation structure (Pulliam 
1988; Hanski 1999). Over time, we may expect the upper 
subpopulations to lose genetic variation and that genetic 
variation is higher at the receiving end of the asymmetric 
gene flow (Hanfling and Weetman 2006). The smablank 
is threatened by habitat modifications due to hydropower 
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development, causing severely reduced water flow and 
habitat fragmentation due to dams and weirs (Thorstad 
et al. 2009). 

On this background, the following hypotheses were 
tested: 

HI: The smablank population exhibits a fine-scale popu- 
lation structure with significant F ST among subpopulations. 

H2: Genetic variation is lower in the upper than in the 
lower region of the distribution area due to expected 
asymmetric gene flow. 

H3: Over the last six decades, which is a period of 
severe habitat reduction and modification, smablank has 
experienced a substantial reduction in genetic variation. 

Material and Methods 
Study area 

River Namsen runs for 210 km in Nord-Trondelag 
County, central Norway, from Lake Namsvatnet (455 m 



a.s.l.) to the outlet in the sea at the town of Namsos 
(64.46°N, 11.51°E) (Fig. 2). The catchment area is 
6,265 km 2 , and mean annual discharge at the outlet to 
the sea is 290 m 3 sec -1 . In the distribution area of 
smablank, the only other major fish species is brown 
trout, whereas three-spined stickleback (Gasterosteus 
aculeatus L.) occurs in the lower areas (A and B, Fig. 2; 
Berg 1984). European minnow (Phoxinus phoxinus L.) has 
been introduced to the upper parts of the river system 
and is spreading downstream, but has not yet reached the 
smablank areas (own unpublished data). 

Sampling of landlocked salmon (smablank) 

Male smablank mature sexually between 2 and 4 years, 
and females between 3 and 5 years. The smallest mature 
males observed were 120 mm long (total length), while 
the smallest mature females were 175 mm (Berg and Gau- 
sen 1988; Thorstad et al. 2009; Norum 2010). The largest 
reported individual was 295 mm long. 
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Figure 2. Location of River Namsen and the 
position of sampling localities A-G. Upper 
panel: map with the original distribution area 
of river-resident Atlantic salmon ("smablank") 
in bold lines. The river section with overlapping 
smablank and anadromous Atlantic salmon is 
indicated by shading. Lower panel: schematic 
gradient of River Namsen. Symbols as in upper 
panel. Marine limit is the highest postglacial 
marine shoreline. 



1540 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



0. T. Sandlund ef al. 



River-Resident Salmon 



Fin-clip samples for genetic analyses were collected 
between 2005 and 2008 from fish caught in six areas in the 
main river (Fig. 2, localities A-F) and from the tributary 
River Mellingselva (Fig. 2, locality G) (Table 1). Within 
each area, fish were collected from several locations in an 
effort to decrease relatedness between individuals. In order 
to analyze the temporal stability of the genetic structure, 
genetic analyses were also performed on scale samples col- 
lected between 1955-1956 and 1978-1979 (in total 100 fish) 
from three river sections (Table 1; identified with the sub- 
scripts 55, 56, 78, and 79). Scales from anadromous Atlantic 
salmon, collected in the lower sections of River Namsen in 
1978 (Table 1, locality ANA 78 , cf. river section ANA in 
Fig. 2), were analyzed in order to compare smablank and 
anadromous salmon from the same river. In total, samples 
from 312 individual fish were subject to genetic analyses. 
The reason for analyzing samples of anadromous salmon 
from 1978 was that they were collected before the start of 
salmon aquaculture. Present-day samples may have been 
influenced by introgression with escaped farmed salmon 
from coastal aquaculture (Glover et al. 2013). 

Distribution and habitat 

The habitat area estimates for smablank were based on 
aerial photos from 3 June 2003 to 22 July 2006 (Thorstad 



et al. 2009) and later field classification of habitat quality 
(Norum 2010). Areas were determined using the GlS-pro- 
gram ArcMap 9.3. 

The original distribution of smablank included the 
main river from above the waterfall Nedre Fiskumfoss 
(elevation top 58 m a.s.l.) up to Namskroken (286 m 
a.s.L), a distance of about 85 km (Berg 1953; Thorstad 
et al. 2009; Fig. 2). Construction of hydropower dams 
and fish passes has allowed access for anadromous Atlan- 
tic salmon up to the waterfall Aunfoss (95 m a.s.L), 
reducing the allopatric distribution of smablank by 10 km 
in the main river. We do not know the status of 
smablank in the area where it is sympatric with anadro- 
mous salmon, and whether smablank and anadromous 
salmon hybridize. Today, the hydropower dam at Aunfoss 
separates anadromous salmon and smablank in the main 
river so that smablank are allopatric in 75 km of the 
main river (Fig. 2, localities A-G) between Aunfoss and 
Namskroken. Along this river section, smablank also 
occur in the lower reaches of a number of tributaries 
(Fig. 2). The total water-covered area where allopatric 
smablank have been recorded in River Namsen and tribu- 
taries amounts to 12.5 km 2 , with the major proportion 
(87%) in the main river. 

The river gradient along the smablank distribution area 
(dammed areas excluded) varies between 0.8 and 



Table 1. Localities, sampling years, and number of samples for genetic analyses of anadromous and river-resident Atlantic salmon (smablank) in 
River Namsen and the tributary River Mellingselva. Sampling localities are indicated in Fig. 2 by their letter symbols: A-G and ANA. Recent samples 
(2005-2008) are fin clips, older samples are scales. Geographical positions of sections are given in decimal degrees from downstream to upstream 
end of actual sampling locality. 



Distance from 



Elevation 



Sampling locality 


Position (decimal degrees) 


Year 




# fish 


sea (km) 


(m a.s.l.) 


ANA 78 


River Namsen, anadromous section 


64,46748°N 11,54462°E- 
64,541 88°N 12,45527°E 


1978 




35 


0-83 


0 


A 


Dam Aunfoss to dam Asmulfoss 


64,61 126°N 12,57973°E- 
64,67698°N 12,67656°E 


2005- 


-08 


14 


84-93 


94 


A 79 




As above 


1979 




26 






B 


Dam Asmulfoss to Trongfoss 


64,67698°N 12,67656°E- 
64,75004°N 12,84708°E 


2005- 


-08 


16 


94-107 


101 


C 


Trongfoss to weir Kjelmyrfoss 


64,75004°N 12,84708°E- 
64,81 540°N 12,96043°E 


2008 




53 


1 08-1 1 6 


144 


D 


Weir Kjelmyrfoss to below weir Bjorhusdal 


64,81 540°N 12,96043°E- 
64,89635°N 13,07082°E 


2005- 


-08 


17 


117-129 


158 


E 


Weir Namsskogan to below Bj0rnstadfoss 


64,931 52°N 13,16340°E- 
65,01 149°N 13,26039°E 


2005- 


-08 


26 


137-147 


210 


E/8 




65,01 149°N 13,26039°E 


1978 




30 










65,01 149°N 13,26039°E 


1955 




15 






F 


Above Bjornstadfoss to outlet River Mellingselva 


65,01555°N 13,25750°E- 
65,05575°N 13,31968°E 


2005- 


-08 


28 


148-153 


221 


G 


River Mellingselva 


65,05575°N 13,31968°E- 
65,07937°N 13,28596°E 


2008 




23 


154 


246 


G 7 8 




As above 


1978 




19 






G 56 




As above 


1956 




10 
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12 m km -1 (Fig. 2). There are several waterfalls, but only 
Trongfoss (cf. Fig. 2), with a vertical fall of more than 

10 m, is considered to be a complete barrier to upstream 
movement. The smaller waterfall Bjornstadfoss (cf. Fig. 2) 
is probably passable for upstream movement at certain 
water discharges. The construction of Asmulfoss dam and 
hydropower station (around 1970) created a permanent 
barrier where only a number of small rapids occurred 
before. Since 1970, weirs have been constructed at several 
sites, acting either as permanent barriers or obstacles at 
certain water discharges. Both dams and weirs cause hy- 
dromorphological changes promoting lentic instead of 
lotic conditions. In total, 38 km of the 75-km-long river 
stretch (51%) available for smablank in the main river 
are presently affected by the change to a more lentic habi- 
tat. In addition to habitat degradation, water has been 
diverted from the main river for hydro-electric power 
production. The effect of this water diversion is that the 
size of the remaining catchment area upstream of Bjorns- 
tadfoss (147 km from the sea; Fig. 1) is presently only 1/3 
of the natural situation. As a consequence, annual water 
discharge and water velocities are reduced. 

Genetic methods 

The fish material was assayed for genetic variation at 
eight microsatellite loci, and a subsample of 88 individu- 
als of smablank from three sampling localities was also 
assayed for genetic variation at 4414 SNP loci (Table 2 
and Fig. 2). The SNP data were available from a previous 
study (Bourret et al. 2013). Although these data concern 
fewer localities and specimens, they were included here to 
add extra statistical power to the analyses. DNA was 
extracted from ethanol preserved fin clips and scales using 
the E.Z.N.A™ tissue DNA kit (E.Z.N.A.® Omega Bio-Tek 
Inc, Norcross, GA). 

PCR was carried out in two different multiplexes con- 
taining the following microsatellite loci: Ssa289, Ssal4 
(McConnell et al. 1995), Ssal71, Ssal97 (O'Reilly et al. 
1996), Ssa408 (Cairney et al. 2000) in multiplex 1, /(20.19 
(Sanchez et al. 1996), Ssosl85 (Slettan et al. 1995), Ssosl438 
(Slettan et al. 1996) in multiplex 2. The multiplex reactions 
were carried out in a total volume of 10 /.(L, containing 

11 /imol/L and 6.5 /miol/L of total primers, in multiplex 1 
and 2, respectively, but with different concentrations for 
each pair, 1 mmol/L of total dNTP, IX reaction buffer, 
2.25 mmol/L MgCl 2 , and 0.75 units of Thermostart taq 
polymerase (Thermo Scientific Inc., Waltham, MA). The 
following PCR program was run on a Quattro Cycler 
(VWR): denaturation for 15 min at 95°C: six cycles touch- 
down PCR of denaturation at 94°C, annealing temperature 
from 58°C to 52 °C, and extension at 72°C for 60 s. The 
last 24 cycles were run with denaturation at 94°C for 30 s, 



Table 2. Summary statistics of eight microsatellite loci and 4414 SNP 
loci from the anadromous and river-resident (smablank) Atlantic sal- 
mon in River Namsen. W is sample size, H e is expected heterozygosity, 
H 0 is observed heterozygosity, #A is observed numbers of different 
alleles, A R is allelic richness based on five diploid individuals, P H w is 
probability of conformance to Hardy-Weinberg expectation, %P is 
proportion of polymorphic SNP loci, A E is average effective number of 
alleles at SNP loci. Sampling localities are indicated in Fig. 2. 



Microsatellite loci 



Sampling locality 


W 


He 




#A 


An 


Ph— w 


ANA 78 


35 


0.722 


0.717 


9.6 


4.9 


0.2802 


A 


14 


0.359 


0.304 


2.9 


2.3 


0.3436 


B 


16 


0.346 


0.344 


2.9 


2.2 


0.9204 


C 


53 


0.377 


0.346 


3.5 


2.4 


0.0042 


D 


17 


0.380 


0.412 


3.5 


2.5 


0.5557 


E 


26 


0.370 


0.418 


3.8 


2.4 


0.1325 


F 


28 


0.361 


0.375 


3.9 


2.6 


0.2479 


G 


23 


0.320 


0.315 


3.4 


2.3 


0.0016 


E55 


15 


0.344 


0.457 


2.8 


2.4 


0.0038 


E/8 


30 


0.355 


0.366 


3.3 


2.5 


0.9849 


G56 


10 


0.305 


0.303 


2.3 


2.1 


0.9544 


G/8 


19 


0.360 


0.322 


3.1 


2.5 


0.1795 


A79 


26 


0.321 


0.298 


3.0 


2.2 


0.2056 




SNP loci 










Sampling locality 


W 


He 


Ho 


%P 


A E 


Ph-w 


ANA 78 


48 


0.345 


0.342 


96.7 


1.591 


-0 


B 


16 


0.126 


0.127 


40.8 


1.213 


~1 


C 


11 


0.120 


0.127 


35.8 


1.204 


-1 


E 


13 


0.128 


0.130 


37.5 


1.217 


~1 



Significant P-values are in bold. 



annealing at 52°C for 30 s, and extension at 72°C for 60 s; 
a final step of extension at 72°C for 10 min. Fragments 
from each multiplex were separated and visualized sepa- 
rately on an ABI 3130x1 DNA analyser (Applied Biosys- 
tems) and sizing using GENEMAPPER ver. 3.7 (Applied 
Biosystems, Foster City, CA). 

SNP genotype data were obtained from two different 
projects. The SNP genotype data for the anadromous 
Atlantic salmon from River Namsen were described by 
Karlsson et al. (2011), using the 7K Atlantic salmon SNP 
chip (CIGENE, As, Norway). The SNP chip data for the 
smablank salmon from River Namsen watercourse were 
described by Bourret et al. (2013), also using the 7K 
Atlantic salmon SNP chip (CIGENE). The two datasets 
were combined, and a common set of 4414 SNPs was 
used in this study. 

Test of conformity to Hardy-Weinberg equilibrium, 
estimates of observed and expected heterozygosity, homo- 
geneity test of allele frequencies, and F ST estimates 
according to Weir and Cockerham (1984) were carried 
out in GENEPOP v.4 (Raymond and Rousset 1995), for 
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microsatellite data and SNP data. Number of alleles and 
number of alleles independent of sample size (allelic rich- 
ness) were estimated for the microsatellite data using 
Fstat v. 2.9.3 (Goudet 2001). Proportion of polymorphic 
loci and effective number of alleles were estimated for the 
SNP data using GENALEX (Peakall and Smouse 2006). 

Possible differences in allelic richness, expected hetero- 
zygosity, and effective number of alleles between sampling 
localities were tested using Wilcoxon signed-ranked test 
as implemented in SPSS Statistics 18 (http://www.spss. 
com/). Possible differences in proportion of polymorphic 
loci (SNPs) were tested by a random resampling 
approach, using POP TOOLS (an add-in program in 
excel available at: http://www.poptools.org/). Random 
sampling was performed 1,000 times, and for each sam- 
pling, the proportion of polymorphic loci was estimated. 
The average proportion of polymorphic loci from each 
random sampling and their upper (0.975) and lower 
(0.025) percentiles were recorded and compared with the 
observed proportion of polymorphic loci in the sampling 
locality for which there were as many individuals as were 
resampled in the locality with a larger sample size. 

Population structure in the smablank was explored from 
pairwise estimates of P S t> analysis of molecular variance 
(AMOVA) implemented in ARLEQUIN ver. 3.5.1.2 
(Schneider et al. 2000), analysis of individual probability 
of belonging to different number of populations without a 
priori information of sampling locality using STRUCTURE 
ver. 2.3.1 (Pritchard et al. 2000), and from pairwise esti- 
mates of genetic distance according to Nei (1972, 1978) 
using Phylip ver. 3.69 (Felsenstein 2009). In STRUCTURE, 
individuals were probabilistically assigned to a predefined 
number of populations, ranging from one (K = 1) to 
twelve (K = 12), without a priori information of sampling 
locality while applying the admixture model. Each run was 
repeated three times with 50,000 repetitions as burn-in 
and 100,000 repetitions after burn-in. The STRUCTURE 
HARVESTER program (Earl and von Holdt 2012) was 
used for estimating the most likely number of populations 
using the Delta_K method (Evanno et al. 2005) and the Ln 
probability of the data estimates (Pritchard et al. 2000). 

Putative first-generation migrants were detected from 
using STRUCTURE ver. 2.3.1 (Pritchard et al. 2000) and 
from individual self-assignment using the Bayesian 
method (Rannala and Mountain 1997) implemented in 
GeneClass2 (Piry et al. 2004). For the latter, the individu- 
als to be assigned were all included in the reference popu- 
lations to which the populations were assigned (self- 
assignment). Genetic self- assignment was also used for 
detecting possible signatures of genetic introgression from 
anadromous salmon in the smablank samples. 

Effective population size was estimated using the tem- 
poral method by Jorde and Ryman (2007), with the Fs 



estimator implemented in the TempoFS software (Jorde 
and Ryman 2007) available at http://www.zoologi.su.se/ 
-ryman/ and the maximum-likelihood method by Wang 
and Whitlock (2003), implemented in the MNe 1.0 
program available at http://www.zsl.org/science/research- 
projects/software/mlne,1151,AR.html. With the latter pro- 
gram, migration rates were estimated using all available 
samples as sources for immigrants. Estimates of effective 
population size using the temporal methods require a 
minimum of two temporal samples from the same popu- 
lation. In this study, temporal samples were available 
from three sampling localities. For the remaining samples 
(and the temporal samples), effective population size was 
estimated using the linkage disequilibrium method (LD) 
by Hill (1981) with the bias correction for sample sizes 
(Waples 2006), as implemented in the LDNe program by 
Waples and Do (2008). The demographic history of the 
populations was examined using the M-test by Garza and 
Williamson (2001). The M-value is the ratio of the num- 
ber of alleles and the allele size range at microsatellite 
loci. A population that has recently experienced a reduc- 
tion in effective population size is expected to show a 
lower M-value, because the number of alleles is expected 
to decline faster than the range in allele size, under a 
step-wise mutation model (smm) (Garza and Williamson 
2001). All microsatellite loci appeared to fit the assump- 
tion of a SMM model, except for one locus (Ssal71) with 
alleles separated by both four and two base pairs. To 
determine whether the estimated M ratio represented a 
signature of a recent decline in population size, the con- 
servative critical M-value of 0.68 was applied, as proposed 
by Garza and Williamson (2001). 

Genetic data 

Of the 312 individuals assayed for genetic variation at 
eight microsatellite loci, 190 had a scoring rate of 100%. 
Genotypes were missing at one locus for 70 individuals, 
at two loci for 27 individuals, at three loci for 14 individ- 
uals, at four loci for six individuals, at five loci for three 
individuals, and at six loci for two individuals. Significant 
deviation from Hardy-Weinberg equilibrium was detected 
in three samples: C, G, and E 55 (Table 2). A closer exami- 
nation of deviation from Hardy-Weinberg equilibrium at 
individual loci and samples detected no consistent trends 
in deficits or excess of heterozygotes (Appendix 1). All 
microsatellite loci were used in the further analyses. 

Of the 88 individuals assayed for genetic variation at 
4414 SNP loci from the locations B, C, and E, 79 individu- 
als had a scoring rate larger than 95%, and all individuals 
had a scoring rate equal to or higher than 90%. Significant 
deviation from Hardy-Weinberg equilibrium was detected 
in the samples from the anadromous population but not in 
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the smablank samples (Table 2). In the anadromous popu- 
lation, 3.7% of the informative loci had a probabil- 
ity < 0.05 of being in Hardy-Weinberg equilibrium. In the 
B, C, and E smablank samples, 4.9, 1.6, and 3.3% of the 
informative loci had probabilities < 0.05 of being in 
Hardy-Weinberg equilibrium, respectively. All SNP loci 
were used in the further analyses. 

Results 

Smablank vs. anadromous Atlantic salmon 

Microsatellites 

Expected heterozygosity and allelic richness within 
smablank sampling locations were much lower (P < 0.012 
for each pairwise comparison, Wilcoxon signed-rank test) 
than the anadromous population (Table 2). Expected het- 
erozygosity in the smablank samples ranged from 0.31 to 
0.38, while expected heterozygosity in the anadromous 
population was 0.72. Allelic richness ranged from 2.1 to 2.6 
in the smablank samples and was 4.9 in the anadromous 
population. Two loci (Ssal4 and |(20.19) were monomor- 
phic in all smablank samples, but polymorphic in the anad- 
romous population, with two and three alleles, respectively 
(Appendix 2). Pairwise estimates of genetic differentiation 
(F sx ) between the anadromous population and the 
smablank samples ranged from 0.22 to 0.26 and were sig- 
nificantly larger than zero (P « 0). There were no signa- 
tures of genetic introgression from anadromous salmon in 
the smablank samples as none of the smablank specimens 
assigned with a relative assignment score higher than 
0.02% to the anadromous sample, whereas all specimens 
had an assignment score -100% to the smablank samples. 

SNPs 

In agreement with the microsatellite data, the anadro- 
mous population had about twice as much genetic varia- 
tion as the smablank samples based on the SNP data 
(Table 2). The expected heterozygosity ranged from 0.12 
to 0.13 in the three smablank sampling localities (B, C, E) 
and was 0.34 in the anadromous population (P « 0 for 
each pairwise comparison, Wilcoxon signed-rank test). 
The anadromous population was polymorphic for 96.6% 
of the loci, while the smablank samples were polymorphic 
for only 35.8-40.8% of the loci. Random sampling (100 
times) of 13 individuals from the anadromous population 
resulted in a mean proportion of polymorphic loci of 
0.91 with a 95% confidence interval of 0.88-0.92, which 
is thus highly significant. 

Pairwise F ST estimates between the anadromous popula- 
tion and the B, C, and E smablank samples were 0.28, 0.27, 



and 0.27, respectively (P « 0). In agreement with the 
results from microsatellite loci, no signatures of genetic 
introgression from anadromous salmon were found in the 
smablank samples, as all smablank specimens had a relative 
assignment score of zero to the anadromous sample, but 
100% to the smablank samples. 

Genetic structure of the smablank 

Microsatellites 

Average expected heterozygosity (range: 0.305-0.377) and 
allelic richness (range: 2.1-2.6) did not differ significantly 
between pairs of smablank sample localities (Wilcoxon 
signed-rank test). 

From pairwise F sr estimates and homogeneity tests 
between pairs of sampling locations (Appendix 2), two 
main groups appeared: Group I comprising the A, B, C, 
and D localities, and Group II comprising the F and G 
localities. The samples from the E locality were most closely 
associated with Group II in 1955-1978 and with Group I in 
2005-2008. The genetic structure was further explored with 
an AMOVA, whereby the E and the E 55 /E 78 samples were 
placed in separate groups or in the two main groups, 
respectively. The largest genetic variation among groups 
(P CT ) was obtained when the E 55 /E 78 samples were placed 
in Group II and the E sample placed in a separate group 
(Per = 0.081, P~0). When the E 55 /E 78 samples and the E 
sample were placed in separate groups, the among-group 
variance component was 0.073 (P~0). When the E 55 /E 78 
samples and the E sample were both placed in Group II, 
the variance among groups was relatively high 
(^ct = 0.080, P = 0.001), but the variance among popula- 
tions within group was larger and significant (P = 0.042). 
From these analyses, there appears to be at least three 
genetic clusters: one in the lower region of the distribution 
area, represented by the samples from A, B, C, and D, one 
in the uppermost region, represented by the samples from 
F and G, and one in between the lower and upper region, 
represented by the E samples. 

The population structure of the smablank samples was 
further explored by individual genetic assignment to a 
predefined number of populations (K) using STRUC- 
TURE. The most likely number of populations was two 
according to the Delta-K value (Delta-K = 25.42) and the 
Ln probability of the data (Mean LnP(K) = -2805.40), 
represented by the A, B, C, and D samples in one cluster, 
the F and G samples in the second cluster, and the E 
samples positioned in between the two clusters (Fig. 3). 

The genetic differences among the smablank sampling 
localities are summarized in a neighbor-joining dendro- 
gram (Fig. 4). The two main clusters, represented by sam- 
ples, form the lower and upper regions of the river are 
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Figure 3. Average proportion of genome membership, from 
genotypes at eight microsatellite loci for individuals of river-resident 
Atlantic salmon (smablank) sampled at eight localities in River 
Namsen, assuming two populations (K = 2) and applying the 
admixture model in STRUCTURE. Sampling localities are A-G (cf. 
Fig. 2). For older samples, sampling years are given as subscript. 
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visible, and the E sample is positioned in between these 
two clusters. The sample collected at B]0rnstad in 1955 
(E 55 ), clustered with the samples from the upper region, 
while the samples collected in 1978 (E 78 ) clustered some- 
what in between the two main clusters. 

SNPs 

There were no significant differences in the effective num- 
ber of alleles, expected heterozygosity (Wilcoxon signed- 
rank test) or in the proportion of polymorphic loci (ran- 
dom resampling) between pairs of sampling localities 
(Table 2). 

The genetic differences indicated by microsatellites, 
between the E sample and the samples from A, B, C, and 
D, were supported by the SNP data with an F ST of 0.046 
(P~0) between the E and B samples, and 0.053 (P~0) 
between the E and C samples. In agreement with the mi- 
crosatellite data, the SNP analysis showed no significant 
difference in allele frequencies between B and C 
(F ST = 0.027, P = 0.99), and the most likely number of 
populations detected by STRUCTURE was two according 
to the Delta-K value (Delta-K = 166.52) but three accord- 
ing to the Ln probability of the data (Mean LnP 
(K) = —68850). From a closer examination of the results 
from the STRUCTURE analyses, assuming three popula- 
tions, (Fig. 5) one population was represented by the 
individuals from E, one population was represented by 
individuals from B and C, and one population was repre- 
sented by only two individuals from B. These two indi- 
viduals had a probability of 0.998 of belonging to the 
third population, which otherwise was not represented in 
our samples. Another individual from B had a probability 
of 0.995 of belonging to the population represented by 
the E samples. The STRUCTURE results were almost 
identical when assuming two instead of three populations. 



Figure 4. Neighbor-joining dendrogram generated from pairwise 
estimates of Nei's genetic distance from samples of river-resident 
Atlantic salmon (smablank) in River Namsen watercourse from eight 
sampling localities, including temporal samples, using eight 
microsatellite loci. Sampling localities are shown in Table 1 and Fig. 2. 




Figure 5. Proportion of genome membership for each individual of 
river-resident Atlantic salmon (smablank) sampled at three localities 
(B, C, and E, separated by blank bars; cf. Fig. 2), in River Namsen 
genotyped at 4414 SNP loci, assuming three populations (K = 3), 
represented by three different colors. 

The exception was that when assuming two populations, 
the two individuals from location B representing the third 
population (K = 3) were assigned to the same popula- 
tions as the individuals from the downstream locations B 
and C. From self-assignment using GeneClass, the one 
individual from B was assigned to the E samples (relative 
log-likelihood score = 100%), in agreement with STRUC- 
TURE, while the other two individuals from the B locality 
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with unknown origin from the STRUCTURE analysis 
were assigned to the B locality from where they were 
sampled (relative log-likelihood score = 100%). These 
observations suggest that there may be three (18.8%) 
first-generation migrants among the individuals sampled 
at the B locality, two for which the populations are not 
included in the SNP dataset, and one with a genetic sig- 
nature matching the E samples. 

Temporal genetic variation in smablank 

Historical samples from the A, E, and G sampling locali- 
ties were analyzed for genetic variation at the eight micro- 
satellite loci (Appendix 1). Average expected 
heterozygosity and allelic richness (Table 2) did not differ 
significantly between historical samples and contemporary 
samples within sampling localities (Wilcoxon signed-rank 
test). No significant temporal instability in allele frequen- 
cies (Appendix 2) was observed between the samples from 
1979 and 2005-2008 at the A locality (F ST = -0.01, 
P = 0.94), nor between samples from 1956, 1978, and 
2008 at the G locality (G 56 vs. G: F sr = 0.032, P = 0.114; 
G vs. G 78 : F ST = 0.011, P = 0.244; G 78 vs. G 56 : 
F S t = 0.034, P = 0.059). However, between the E samples 
from 1955 and 1978, a significant temporal instability in 
allele frequencies was observed (F ST = 0.04, P = 0.002). 

The impression of a general trend of decreasing genetic 
differentiation among localities over time emerging from 
Fig. 2 was specifically tested for pairs of localities where 
temporal samples were available, that is, between A and E, 
A and G, and E and G (cf. Fig. 2). The genetic difference 
between the A and E was highly significant (P = 0.000) in 
1979/1978 and significant (P = 0.020) in 2005/2008, and 
the estimated F sr was higher in 1978/1979 (F ST = 0.084) 
compared with the 2005/2008 samples (F ST = 0.018). The 
genetic difference between the A and G localities was highly 
significant (P = 0.000) in 1979/1978 and also in 2005/2008, 
but the F S t estimate was somewhat higher in 1978/1979 
(Fgr = 0.112) than in 2005/2008 (F sr = 0.077). The genetic 
difference between the E and G localities was higher and 
significant in the 1955/1956 samples (F ST = 0.051, 
P = 0.001), compared with the nonsignificant 1978 samples 
(F ST = 0.019, P = 0.076). The 2005/2008 samples from the 
E and G localities were significantly different (P = 0.021), 
but with an F ST estimate (F ST = 0.02) not significantly dif- 
ferent from 1978. Jack-knifing over five polymorphic loci 
demonstrated that the latter test had low power (i.e., large 
estimated standard errors). 

Effective population size and migration rate 

Based on the temporal methods, effective smablank popu- 
lation size (N e ) ranged from 29 (E locality) to 302 (G 



Table 3. Estimates of effective population size (Ne) using temporal 
methods (TempoFS and MINe) and the linkage disequilibrium method 
(LDNe) for samples of anadromous and river-resident (smablank) 
Atlantic salmon in River Namsen. Estimates of migration rates and Ne 
were jointly obtained from the temporal method by Wang and Whit- 
lock (2003). One dataset with eight microsatellite loci and one dataset 
with 4414 SNP loci were used as indicated in the table (Marker type). 

Temporal method - Tempo FS 



Pop 


Ne (95% CI) 




Marker type 


E55-E78 


29 (16-168) 




Msat 


G56-G78 


261 (28-oo) 




Msat 


G78-G 


88 (27-oo) 




Msat 


G 56 -G 


126 (35-00) 




Msat 


A79-A 


cc (178-oo) 




Msat 


Temporal method - MINe 






Pop 


Ne (95% CI) 


m /o \-\) 


Marker type 


E55-E78 


36 (14-176) 


0.12 (0.042-0.35) 


Msat 


G56-G78 


30 (9-oo) 


0.99 (0.09-oo) 


Msat 


G78-G 


302 (54-oo) 


0.006 (0.001-0.045) 


Msat 


G56-G 


79 (28-373) 


0.038 (0.0062-0.1) 


Msat 


A79-A 


213 (77-oo) 


0.17 (00-00) 


Msat 


Moment estimate - LDNe 






Pop 


Ne (95% CI) 




Marker type 


ANA 78 


00 (141-oo) 




Msat 


A, B, C, D 


266 (54-oo) 




Msat 


E 


00 (106-oo) 




Msat 


F, G 


58 (22-3216) 




Msat 


ANA 78 


268 (263-274) 




SNPs 


B, C 


28 (27-28) 




SNPs 


E 


619 (406-1295) 




SNPs 



locality), with large confidence intervals for some esti- 
mates (Table 3). The estimate of migration rates was 
0.038 in the G locality, while increasing further down- 
stream in the E locality (0.12) and the A locality (0.17), 
with large confidence intervals (Table 3). The linkage dis- 
equilibrium method gave poor results based on microsat- 
ellite data, with negative estimates and very large or 
infinite confidence intervals. Based on SNP data, N e was 
268 in the anadromous sample, while N e for smablank 
was 28 at the B and C localities combined, and 619 for 
the E locality (Table 3). The larger N e estimate for the E 
locality compared with the anadromous samples might 
reflect a possible bias (Waples and Do 2010) and repre- 
sent the effective population size of the metapopulation 
including localities E, F, and G (cf. Gomez-Uchida et al. 
2013). All smablank samples had significantly lower M- 
values than the conservative critical M-value (0.68) 
proposed by Garza and Williamson (2001). The average 
M-value for the anadromous population was 0.87 using 
seven informative microsatellite loci (range: 0.55-1.0). For 
the smablank samples, there were five informative loci. 
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The average M-value for the pooled contemporary sam- 
ples from A, B, C, and D was 0.54 (range: 0.3-0.70). For 
the E sample, the average M-value was 0.49 (range: 0.22- 
0.75). The average M-value for the pooled contemporary 
samples from F and G was 0.55 (range: 0.33-0.8). M-val- 
ues for the historical samples did not appear to differ 
from the contemporary samples, with M = 0.53 for the 
Ess/Eye samples and 0.42 for the A 79 sample. The low M- 
values in smablank should not immediately be interpreted 
as a sign of a recent bottleneck because the M-values 
might be biased from a few outlier alleles (especially in 
the Ssosl438 and Ssosl85 loci) and because only five loci 
were examined. Furthermore, the number and origin of 
founders, the long-term effective population size (N e ), 
and possible fluctuations in N e are not known, which 
makes it difficult to deduce whether the M-values are due 
to recent or past demographic events (Garza and Wil- 
liamson 2001). 

Discussion 

The heterozygosity and allelic richness of the river-resident 
smablank was about 50% lower than that of the anadro- 
mous Atlantic salmon in the lower river sections. We 
found no signature of genetic introgression from releases 
of offspring of anadromous salmon in the smablank sam- 
ples. These findings confirm earlier observations in 
enzyme-coding genes for smablank (Vuorinen and Berg 
1989). Smablank have been isolated from other conspecific 
populations for approximately 9,500 years (Berg 1985), 
corresponding to some 2000 generations. It is reasonable 
to assume that a major reduction in genetic variation was 
caused by founder effects and genetic drift during the first 
generations after isolation. The processes when smablank 
became isolated from anadromous Atlantic salmon several 
thousand years ago are not known, but the low frequency 
of mature resident females in anadromous populations 
(Hindar and Nordland 1989) indicates that the founding 
event may have involved only a few resident females. 
Genetic drift and other stochastic processes in restricted 
populations may have caused the low genetic diversity 
indicated by the low heterozygosity in smablank as well as 
among other landlocked populations compared with anad- 
romous populations (Vuorinen and Berg 1989; Ozerov 
et al. 2012; Bourret et al. 2013; Perrier et al. 2013). The 
level of genetic variation in the smablank (island) popula- 
tion relative to the anadromous salmon (mainland) popu- 
lation is similar to island/mainland pairs in other taxa that 
have been separated for a similar time period (e.g., Frank- 
ham 1997; Hinten et al. 2003; Whiteley et al. 2010). 

A large-scale genetic comparison based on SNPs of 
38 Atlantic salmon populations in Europe and North 
America revealed that nonanadromous populations repre- 



sented distinct outliers in all regions (Bourret et al. 2013). 
Among anadromous populations, differentiation between 
European and North American populations accounted for 
41% of variation, whereas within Europe, three major 
groups (Atlantic, Baltic, and Barents-White Sea) 
accounted for 8% of the variation. Smablank clustered dif- 
ferently from anadromous salmon of European origin and, 
genetically, was no more closely related to Norwegian than 
to other European anadromous populations (Fig. 2 and 
Table SI in Bourret et al. 2013). Based on the 4414 SNP 
loci used in the present study, estimated E S t between the 
smablank samples and the anadromous population in 
River Namsen ranged from 0.27 to 0.28, while estimated 
f ST among 13 Norwegian anadromous populations using 
the same SNP loci was 0.038 (Karlsson et al. 2011). Based 
on microsatellite loci, the genetic distance between 
smablank and anadromous Atlantic salmon in River Nam- 
sen was estimated at F ST = 0.24, which is similar to the 
difference between anadromous Atlantic salmon in Europe 
and North America (Stahl 1987; King et al. 2007). 

Our results indicate that smablank is subdivided into 
three genetically different populations. One population 
occupies the habitat from Aunfoss up to the weir at Kjel- 
myrfoss (section A-D, a river section of more than 
30 km). The second population inhabits the river section 
from the weir at Namsskogan up to the Bjornstadfoss 
waterfall (section E, approximately 10 km). The third 
population occupies the main river from the Bj0rnstad- 
foss waterfall up to and including the tributary River Mel- 
lingselva (section F and G, a little more than 10 km). The 
subdivision of smablank in several populations based on 
samples from 2005-2008 was also indicated in the sam- 
ples from the 1950s and 1970s. Similarly, enzyme electro- 
phoresis performed in the 1980s indicated a subdivision 
of the smablank in two different genetic clusters corre- 
sponding to the upper and lower regions of the distribu- 
tion (Vuorinen and Berg 1989). This subdivision is 
maintained within a limited habitat area of 12.5 km 2 , and 
it is likely determined by obstacles to upstream move- 
ment. The genetic differentiation is to some extent 
reflected in life-history differentiation. Body size at matu- 
rity is smaller in fish sampled in the upper area (localities 
F and G), compared with other localities (Berg and Gau- 
sen 1988; Thorstad et al. 2009). The morphology of sub- 
populations has not been investigated. 

Effective population size estimates indicate N e values at 
a few thousand smablank in total. Assuming that the divi- 
sion in subpopulations is a persistent pattern, each sub- 
population is at N e values of a few hundred. As expected, 
only downstream migrants were detected by genetic 
analyses, suggesting asymmetric gene flow. This has likely 
occurred ever since smablank became isolated as a land- 
locked population and caused the establishment of a 
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source-sink metapopulation structure. This is expected in 
fast-flowing rivers where it is easier to migrate down- 
stream than upstream (Kawecki and Holt 2002). 

Our expectation that genetic variation was lower in the 
upper (F, G) than in the lower (A-D) subpopulations 
was not confirmed. Neither had the relative level of 
genetic variation in the upper and lower subpopulations 
changed over the last six decades, since the 1950s. The 
genetic variation in the source population seems not to 
be declining, and it is not lower than the sink popula- 
tions, in spite of the asymmetric downstream gene flow. 
This may indicate an effective source population size (N e ) 
large enough to prevent loss of genetic variation in a 
short period of time (since the 1950s), and larger than 
the effective sink population size. It should be noted that 
the area represented by the G samples (the River Melling- 
selva tributary) is in a more or less pristine condition in 
terms of human encroachment. The time since the 1950s, 
corresponding to 8-12 smablank generations, may also 
have been too short for detecting loss of genetic variation. 
Moreover, we cannot exclude the possibility that 
smablank had already been impacted by hydropower 
development when the first samples were collected in 
1955-1956 and 1978-1979. The first hydropower develop- 
ment to influence a major part of River Namsen was the 
damming of the Namsvatn reservoir (the source of River 
Namsen) in 1952. Both during the construction and the 
operational phase, this would have impacted on a major 
part of the smablank habitat. The weirs in the smablank 
section of River Namsen were constructed in 1965/1966, 
1978, and 1998 (Thorstad et al. 2009). 

The subpopulations of smablank are genetically diver- 
gent, but our data indicate that they over the last decades 
have become more similar. This might be caused by an 
increase in asymmetrical gene flow where downstream 
movement has become more frequent relative to 
upstream movement. Whether a downstream sink popu- 
lation depends on the upstream source population for its 
persistence depends, i.a., on the effective population size 
of the sink population and different habitat qualities 
between the source and the sink population (Kawecki and 
Holt 2002). The construction of weirs has likely impacted 
more on upstream than on downstream movement rates, 
thereby increasing the asymmetry of the gene flow. The 
modifications of the river habitat have been more exten- 
sive in the lower part of the distribution area. 

The genetic division in subpopulations is probably an 
advantage to the conservation of genetic variation in 
smablank, because a number of subpopulations with lim- 
ited gene flow between them may conserve more genetic 
variation than a similar number of individuals in one 
continuous population (Hanski 1999; Tufto and Hindar 
2003; Schindler et al. 2010). This is true as long as the 



effective population size is large enough, and limited 
movement and gene flow are possible between subpopula- 
tions to prevent them from extinction. Human activities 
leading to fragmentation of the smablank habitat may 
have reduced or even terminated gene flow among sub- 
populations. Such fragmentation is expected to cause 
increased loss of genetic variation, increased inbreeding, 
and increased risk of extinction (Hanski 1999; Morita and 
Yamamoto 2002; Hanfling and Weetman 2006). To main- 
tain a healthy metapopulation of smablank, measures 
should be taken to maintain a natural pattern of two-way 
movement among subpopulations. Further construction of 
weirs and other movement barriers is not recommended. 

More than half of the distribution area of smablank in 
the main river has over the last decades been converted 
from lotic to lentic habitat due to the construction of 
dams and weirs. The weirs are low dams constructed in 
order to increase water-covered areas at reduced water 
flows, largely for esthetic purposes and to improve access 
to trout fishing. The dams and weirs contribute to 
remove the stony habitats with fast current used by 
smablank (Norum 2010), which seem to have similar 
habitat requirements as juvenile anadromous Atlantic sal- 
mon (cf. Heggenes et al. 1999). The fact that about 60% 
of the annual water volume has been diverted from the 
upper part of the main river contributes to a reduced 
annual flow and reduced water velocities. This impacts 
approximately 53 km, or more than 70%, of the section 
of the main river inhabited by allopatric smablank. In 
addition to the increased area with lentic characteristics, 
the regulated and stabilized water flow of the main river 
probably contributes to sediment packing of the stony 
substratum and reduced shelter availability (Finstad et al. 
2007). Hence, tributaries may presently become more 
important as habitat for smablank, as they are less 
impacted by anthropogenic activities. However, the area 
available to smablank in the tributaries is restricted to 
only 2.1 km 2 (Norum 2010). More detailed sampling is 
required to provide a better understanding of the fine- 
scale population structure and to identify the most 
important areas for smablank. 

Many of the European landlocked salmon populations 
have been subject to heavy pressures due to anthropo- 
genic activities (Berg 1985; Barlaup et al. 2005; Ozerov 
et al. 2012). Although this is also the case for smablank, 
our hypothesis that these anthropogenic changes were 
reflected in a loss of genetic variation (heterozygosity and 
allelic richness) over the last five to seven decades was 
not supported. 

Four features are usually involved in the concept of rar- 
ity: number of individuals and populations, geographical 
distribution, ecological divergence (habitat specificity), and 
genetic divergence (Groom et al. 2006). Contrary to the 
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river-resident Atlantic salmon in North America (Scott 
and Crossman 1973), which occur in several river basins 
(Webb et al. 2007), smablank is the only one in Europe. 
The river-resident populations in North America belong to 
another genetic subgroup than smablank (Lubieniecki 
et al. 2010; Bourret et al. 2013). Smablank clearly represent 
a unique genetic cluster, in line with the concept of Evolu- 
tionary Significant Unit (ESU), because it is matched by no 
other Atlantic salmon population, neither in biologic nor 
genetic characteristics (Waples 1991; Ford 2004). 

In conclusion, our results showed a relatively low genetic 
variation within smablank compared with anadromous 
Atlantic salmon, but still we found a clear genetic subpop- 
ulation structure. Hence, our first hypothesis was partly 
supported. The observed genetic variation and the com- 
plete isolation make smablank an endemic island popula- 
tion. Because the extinction risk is higher in island than 
mainland populations (Frankham 1997), special care is 
needed in the management to conserve such populations. 

A downstream asymmetric gene flow was recorded 
within smablank, but the hypothesis predicting a larger 
genetic variation in the lower compared with the upper 
subpopulations was not supported. This may indicate that 
populations in the pristine habitat in the upper section of 
the smablank area act as a source population, while the 
populations in the lower section are dependent on the 
upper populations to maintain their genetic variation. 

Our third hypothesis was not supported, as we found 
no reduction in genetic variation in smablank over the 
last 50-60 years (12-15 generations). However, we 
observed a reduced differentiation among the subpopula- 
tions within the distribution area. 

A major reason for the apparent genetic stability of 
smablank is likely that its metapopulation structure is still 
relatively intact. To reduce the risk of extinction and 
ensure the viability of this endemic form of Atlantic sal- 
mon, measures causing further fragmentation of the pop- 
ulation should be avoided. 

Acknowledgments 

We gratefully acknowledge Torveig Balstad and Line Erik- 
sen (NINA) for DNA extraction and microsatellite geno- 
typing. Assistance during collection of the contemporary 
samples was provided by Laila Saksgard, Trygve Hestha- 
gen, and Anders G. Finstad. The study was funded by 
NINA and NTNU, with some contributions from the 
County Governor of Nord-Trondelag and the hydropower 
company Nord-Trandelag Elektrisitetsverk (NTE). 

Data Accessibility 

Microsatellite data information: doi:10.5061/dryad.jk42p 

© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



Conflict of Interest 

None declared. 
References 

Barlaup, B. T., E. Kleiven, H. Christensen, N. B. Kile, B. O. 

Martinsen, and A. Vethe. 2005. Bleka i Byglandsfjorden - 

bestandsstatus og tiltak for 0kt naturlig rekruttering. 

DN-utredning 2005-3.72 s. Direktoratet for 

naturforvaltning, Trondheim, Norway. (In Norwegian). 
Berg, M. 1953. A relict salmon, Salmo salar L., called 

"smablank" from the River Namsen, North-Trondelag. Acta 

Bor. A. Sci. 6:1-17. 
Berg, O. K. 1984. comparison between the distributions of 

land-locked Atlantic salmon Salmo salar L. and three-spines 

stickleback Gasterosteus aculeatus L. in the river Namsen, 

Norway. Fauna norv. Ser. A 5:37-41. 
Berg, O. K. 1985. The formation of non-anadromous 

populations of Atlantic salmon, Salmo salar L., in Europe. J. 

Fish Biol. 27:805-815. 
Berg, O. K., and D. Gausen. 1988. Life history of a riverine, 

resident Atlantic salmon Salmo salar L. Fauna Norvegica Ser 

A 9:63-68. 

Bourret, V., M. P. Kent, C. R. Primmer, A. Vasemagi, S. 
Karlsson, K. Hindar, et al. 2013. SNP-array reveals 
genome-wide patterns of geographical and potential 
adaptive divergence across the natural range of Atlantic 
salmon (Salmo salar). Mol. Ecol. 22:532-551. 

Bronmark, C, and L.-A. Hansson. 2005. The biology of ponds 
and lakes. 2nd ed. Oxford Univ. Press, Oxford. 

Cairney, M., J. B. Taggart, and B. Hoyheim. 2000. 

Characterization of microsatellite and minisatellite loci in 
Atlantic salmon (Salmo salar L.) and cross-species 
amplification in other salmonids. Mol. Ecol. 9:2175-2178. 

Dahl, K. 1928. The "blege" or dwarf-salmon. A landlocked 
salmon from Lake Byglandsfjorden in Setesdal. Skrifter, 
Norske Videnskaps-Akademi, Oslo. 
Matematisk-Naturvitenskapelig Klasse 9, 1927. 28 pp. 

Earl, D. A, and B. M. von Holdt. 2012. STRUCTURE 
HARVESTER: a website and program for visualizing 
STRUCTURE output and implementing the Evanno 
method. Conserv. Genet. Resour. 4:359-361. 

Evanno, G., S. Regnaut, and J. Goudet. 2005. Detecting the 
number of clusters of individuals using the software 
STRUCTURE: a simulation study. Mol. Ecol. 14:2611-2620. 

Felsenstein, J. 2009. PHYLIP (Phylogeny Inference Package) 
version 3.69. Distributed by the author. Department of 
Genetics, University of Washington, Seattle. 

Finstad, A. G., S. Einum, T. Forseth, and O. Ugedal. 2007. Shelter 
availability affects behaviour, size-dependent and mean growth 
of juvenile Atlantic salmon. Freshw. Biol. 52:1710-1718. 

Ford, M. J. 2004. Conservation units and preserving diversity. 
Pp. 338-357 in A. P. Hendry and S. C. Stearns, eds. 



1549 



River-Resident Salmon 



O. T. Sandlund ef al. 



Evolution illuminated. Salmon and their relatives. Oxford 
Univ. Press, New York. 

Frankham, R. 1997. Do island populations have less genetic 
variation than mainland populations? Heredity 78:311-327. 

Garza, J. C, and E. G. Williamson. 2001. Detection of 

reduction in population size using data from microsatellite 
loci. Mol. Ecol. 10:305-318. 

Gibson, R. J., D. D. Williams, C. McGowan, and W. S. 
Davidson. 1996. The ecology of dwarf fluvial Atlantic 
salmon, Scdmo salar L., cohabiting with brook trout, 
Salvelinus fontinalis (Mitchill), in southeastern 
Newfoundland, Canada. Pol. Arch. Hydrobiol. 43:145-166. 

Glover, K. A., C. Pertoldi, F. Besnier, V. Wennevik, M. Kent, 
and 0. Skaala. 2013. Atlantic salmon populations invaded 
by farmed escapees: quantifying genetic introgression with a 
Bayesian approach and SNPs. BMC Genet. 14:74. doi:10. 
1186/1471-2156-14-74 

Gomez-Uchida, D., F. P. Palstra, T. W. Knight, and D. E. 
Ruzzante. 2013. Contemporary effective population and 
metapopulation size (N e and meta-jV e ): comparison among 
three salmonids inhabiting a fragmented system and 
differing in gene flow and its asymmetries. Ecol. Evol. 
3:569-580. 

Goudet, J. 2001. FSTAT, a program to estimate and test gene 

diversities and fixation indices (version 2.9.3). Available at 

http://www.unil.chlizea/softwares/fstat.html. 
Groom, M. J., G. K. Meffe, and C. R. Carroll. 2006. Principles 

of conservation biology, 3rd ed. Sinauer, Sunderland, MA. 
Hanfling, B., and D. Weetman. 2006. Concordant genetic 

estimators of migration reveal anthropogenically enhanced 

source-sink population structure in the river sculpin, Cottus 

gobio. Genetics 173:1487-1501. 
Hanski, I. 1999. Metapopulation ecology. Oxford Univ. Press, 

Oxford. 

Heggenes, J., J. L. Bagliniere, and R. A. Cunjak. 1999. Spatial 

niche variability for young Atlantic salmon (Salmo salar) 

and brown trout (S. trutta) in heterogeneous streams. Ecol. 

Freshw. Fish 8:1-21. doi:10.1111/j.l600-0633.1999.tb00048.x 
Hill, W. G. 1981. Estimation of effective population size from 

data on linkage disequilibrium. Genet. Res. 38:209-216. 
Hindar, K., and J. Nordland. 1989. A female Atlantic salmon, 

Salmo salar L., maturing sexually in the parr stage. J. Fish 

Biol. 35:461-463. 
Hinten, G., F. Harriss, M. Rossetto, and P. R. Braverstock. 

2003. Genetic variation and island biogeography: 

microsatellite and mitochondrial DNA variation in island 

populations of the Australian bush rat, Rattus fuscipes greyii. 

Conserv. Genet. 4:759-778. 
Jonsson, B., and N. Jonsson. 2011. Ecology of Atlantic salmon 

and brown trout. Habitat as a template for life histories. 

Springer, Dordrecht. 
Jorde, P. E., and N. Ryman. 2007. Unbiased estimator for 

genetic drift and effective population size. Genetics 177:927- 

935. 



Karlsson, S., T. Moen, S. Lien, and K. Hindar. 2011. Generic 
genetic differences between farmed and wild Atlantic salmon 
identified from a 7K SNP-chip. Mol. Ecol. Resour. ll(Suppl 
l):247-253. 

Kawecki, T. J., and R. D. Holt. 2002. Evolutionary 
consequences of asymmetric dispersal rates. Am. Nat. 
160:333-347. 

Kazakov, R. V. 1992. Distribution of Atlantic salmon, Salmo salar 
L., in freshwater bodies of Europe. Aquae. Res. 23:461-475. 

King, T. I., E. Verspoor, A. P. Spindle, R. Gross, R. B. Phillips, 
M. I. Koljonen, et al. 2007. Biodiversity and population 
structure. Pp. 117-166 in E. Verspoor, L. Stradmeyer and 
J. Nielsen, eds. The Atlantic salmon: genetics, conservation, 
and management. Blackwell, London. 

Lubieniecki, K., S. L. Jones, E. A. Davidson, J. Park, B. F. 
Koop, S. Walker, et al. 2010. Comparative genomic analysis 
of Atlantic salmon, Salmo salar, from Europe and North 
America. BMC Genet. 11:105. doi:10.1186/1471-2156-ll-105 

McConnell, S., L. Hamilton, D. Morris, D. Cook, D. Paquet, 
P. Bentzen, et al. 1995. Isolation of salmonid microsatellite 
loci and their application to the population genetics of 
Canadian east coast stocks of Atlantic salmon. Aquaculture 
137:19-30. 

Morita, K., and S. Yamamoto. 2002. Effects of habitat 
fragmentation by damming on the persistence of 
stream-dwelling charr populations. Conserv. Biol. 16:1318- 
1323. 

Nei, M. 1972. Genetic distance between populations. Am. Nat. 
106:283-392. 

Nei, M. 1978. Estimation of average heterozygosity and genetic 
distance from a small number of individuals. Genetics 
76:379-390. 

Nilsson, J., R. Gross, T. Asplund, O. Dove, H. Jansson, J. 
Kelloniemi, et al. 2001. Matrilinear phylogeography 
of Atlantic salmon (Salmo salar L.) in Europe and 
postglacial colonization of the Baltic Sea area. Mol. Ecol. 
10:89-102. 

Norum, I. C. J. 2010. Habitatkrav og habitattilgjengelighet for 
smablank (Salmo salar), relikt laks i 0vre Namsen. 
Cand.scient. [M.Sc. thesis], Norwegian University of Science 
and Technology, Trondheim, 36 pp + attachments. (In 
Norwegian). 

O'Reilly, P. T., L. C. Hamilton, S. K. McConnell, and J. M. 
Wright. 1996. Rapid analysis of genetic variation in Atlantic 
salmon (Salmo salar) by PCR multiplexing of dinucleotide 
and tetranucleotide microsatellites. Can. J. Fish. Aquat. Sci. 
53:2292-2298. 

Ozerov, M. Y., A. J. Veselov, J. Lumme, and C. R. Primmer. 
2010. Genetic structure of freshwater Atlantic salmon 
(Salmo salar L.) populations from the lakes Onega and 
Ladoga of northwest Russia and implications for 
conservation. Conserv. Genet. 11:1711-1724. 

Ozerov, M. Y., A. E. Veselov, J. Lumme, and C. R. Primmer. 
2012. "Riverscape" genetics: river characteristics influence 



1550 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



0. T. Sandlund ef al. 



River-Resident Salmon 



the genetic structure and diversity of anadromous and 
freshwater Atlantic salmon {Salmo salar) populations in 
northwest Russia. Can. J. Fish. Aquat. Sci. 69:1947-1958. 

Peakall, R., and P. E. Smouse. 2006. GENALEX 6: genetic 
analysis in excel. Population genetic software for teaching 
and research. Mol. Ecol. Notes 6:288-295. 

Perrier, C, V. Bourret, M. P. Kent, andL. Bernatchez. 2013. Parallel 
and non-parallel genome-wide divergence among replicate 
population pairs of freshwater and anadromous Atlantic 
salmon. Mol. Ecol. 22:5577-5593. doi: 10. 1 1 1 1/mec. 12500 

Piry, S., A. Alapetite, J. M. Cornuet, D. Paetkau, L. Baudouin, 
and A. Estoup. 2004. GeneClass2: a software for genetic 
assignment and first-generation migrant detection. J. Hered. 
95:536-539. 

Pollux, B. J. A., A. Luteijn, J. M. Van Groenendael, and N. J. 

Ouborg. 2009. Gene flow and genetic structure of the 

aquatic macrophyte Sparganium emersum in a linear 

unidirectional river. Freshw. Biol. 54:64-76. 
Pringle, J. M., A. M. H. Blakeslee, J. E. Byers, and J. Roman. 

2011. Asymmetric dispersal allows an upstream region to 

control population structure throughout a species' range. 

Proc. Natl Acad. Sci. 108:15288-15293. 
Pritchard, J. K., M. Stephens, and P. Donnelly. 2000. Inference 

of population structure using multilocus genotype data. 

Genetics 155:945-959. 
Pulliam, H. R. 1988. Sources, sinks, and population regulation. 

Am. Nat. 132:652-661. 
Rannala, B., and J. L. Mountain. 1997. Detecting immigration 

by using multilocus genotypes. Proc. Natl Acad. Sci. 

94:9197-9201. 

Raymond, M., and F. Rousset. 1995. Genepop (version 2.1): 
Population genetics software for exact tests and 
ecumenicism. J. Hered. 86:248-249. 

Saisa, M., M.-L. Koljonen, R. Gross, J. Nilsson, J. Tahtinen, J. 
Koskiniemi, et al. 2005. Population genetic structure and 
postglacial colonization of Atlantic salmon {Salmo salar) in 
the Baltic Sea area based on microsatellite DNA variation. 
Can. J. Fish. Aquat. Sci. 62:1887-1904. doi:10.1139/F05-094 

Sanchez, J. A., C. Clabby, D. Ramos, G. Blanco, F. Flavin, E. 
Vazquez, et al. 1996. Protein and microsatellite single locus 
variability in Salmo salar L. {Atlantic salmon). Heredity 
77:423-432. 

Schindler, D. E., R. Hilborn, B. Chasco, C. P. Boatright, T. P. 
Quinn, L. A. Rogers, et al. 2010. Population diversity and 
the portfolio effect in an exploited species. Nature 465:609- 
612. doi:10.1038/nature09060 

Schneider, S., D. Roessli, and L. Excoffier. 2000. Arlequin Ver. 
2.0: a software for population genetic data analysis. Genetics 
and Biometry Laboratory, University of Geneva, Geneva. 

Scott, W. B.,andE. J. Crossman. 1973. Freshwater fishes of Canada. 
Fish. Res. Bd. Can.; Bull. 184: 1-966. 

Slettan, A., I. Olsaker, and 0. Lie. 1995. Atlantic salmon, 
Salmo salar, microsatellites at the SSOSL25, SSOSL85, 
SSOSL311, SSOSL417 loci. Anim. Genet. 26:281-282. 



Slettan, A., I. Olsaker, and 0. Lie. 1996. Polymorphic 
Atlantic salmon, Salmo salar, microsatellites at the 
SSOSL438, SSOSL439 and SSOSL444 loci. Anim. Genet. 
27:57-58. 

Stahl, G. 1987. Genetic population structure of Atlantic 

salmon. Pp 121-140 in N. Ryman, F. Utter, eds. Population 
genetics and fishery management. Univ. of Washington 
Press, Seattle, WA. 

Stahl, G, and K. Hindar. 1988. Genetisk struktur hos norsk 
laks: status og perspektiver. Rapport 1-1988, 
Fiskeforskningen, Direktoratet for naturforvaltning, 
Trondheim, 57 pp. (In Norwegian). 

Thorstad, E. B., K. Hindar, O. K. Berg, L. Saksgard, I. C. J. 
Norum, O. T. Sandlund, et al. 2009. Status for 
smablankbestanden i Namsen. NINA Rapport 403:1-95 (In 
Norwegian). 

Tufto, J., and K. Hindar. 2003. Effective size in management 
and conservation of subdivided populations. J. Theor. Biol. 
222:273-281. 

Vuorinen, J., and O. K. Berg. 1989. Genetic divergence of 
anadromous and nonanadromous Atlantic salmon {Salmo 
salar) in the River Namsen, Norway. Can. J. Fish. Aquat. 
Sci. 46:406-409. 

Wang, J., and M. C. Whitlock. 2003. Estimating effective 
population size and migration rates from genetic samples 
over space and time. Genetics 163:429-446. 

Waples, R. S. 1991. Pacific salmon, Oncorhynchus spp., and the 
definition of "species" under the Endangered Species Act. 
Mar. Fisheries Rev. 53:11-22. 

Waples, R. S. 2006. A bias correction for estimates of effective 
population size based on linkage disequilibrium at unlinked 
gene loci. Conserv. Genet. 7:167-184. 

Waples, R. S., and C. Do. 2008. LDNE: a program for 

estimating effective population size in species from data on 
linkage disequilibrium. Mol. Ecol. Res. 8:753-756. 

Waples, R. S., and C. Do. 2010. Linkage disequilibrium 
estimates of contemporary Ne using highly variable 
genetic markers: a largely untapped resource for 
applied conservation and evolution. Evol. Appl. 3:244- 
262. 

Webb, J., E. Verspoor, N. Aubin-Horth, A. Romakkaniemi, 

and P. Amiro. 2007. The Atlantic salmon. Pp. 17-56 in E. 

Verspoor, L. Stradmeyer and J. Nielsen, eds. The Atlantic 

salmon. Blackwell, Oxford. 
Weir, B. S., and C. C. Cockerham. 1984. Estimating F-statistics 

for the analysis of population structure. Evolution 38:1358- 

1370. 

Whiteley, A. R., K. Hastings, J. K. Wenburg, C. A. Frissell, J. 
C. Martin, and F. W. Allendorf. 2010. Genetic variation and 
effective population size in isolated populations of coastal 
cutthroat trout. Conserv. Genet. 11:1929-1943. 

Wootton, R. J. 1998. Ecology of teleost fishes, 2nd ed. Kluwer, 
London. 



2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



1551 



River-Resident Salmon 



O. T. Sandlund ef al. 



Appendix 1. Summary statistics for eight microsatellite loci for one anadromous Atlantic salmon population and eight sampling 
localities of the landlocked Atlantic salmon in the River Namsen watercourse in Norway. 
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Appendix 1. Continued. 



Locus 



Ssa289 



SsOSL438 



SsOSL85 



Ssa14 



Ssa171 



U20.19 



Ssa197 



Ssa408 



G 

N 

#A 

A R 

Ho 

He 

Ph-w 
E55 
N 
#A 
Ar 
Ho 
He 

Ph-w 
E 78 
N 
#A 
A R 
Ho 
He 

Ph-w 
G 56 
N 
#A 
Ar 
Ho 
He 

Ph-w 

G/8 

N 
#A 

Ar 

Ho 
He 

Ph-w 
A79 
N 
#A 

Ar 

Ho 
He 

Ph w 



23 
4 

2.334 
0.261 
0.303 
0.009 

15 
2 

1.719 
0.200 
0.180 
0.667 

30 
3 

2.217 
0.367 
0.313 
0.680 

10 
1 

1.000 
0.000 
0.000 
na 

19 
3 

2.267 
0.368 
0.314 
0.809 

26 
4 

2.379 
0.615 
0.510 
0.704 



23 
2 

1.217 
0.043 
0.043 
0.915 

14 
1 

1.000 
0.000 
0.000 
na 

26 
2 

1.481 
0.115 
0.109 
0.755 

10 
1 

1.000 
0.000 
0.000 
na 

19 
1 

1.000 
0.000 
0.000 
na 

23 
3 

1.609 
0.130 
0.124 
0.990 



23 
2 

1.792 
0.174 
0.227 
0.263 

13 
2 

1.983 
0.538 
0.393 
0.184 

16 
2 

1.918 
0.250 
0.305 
0.473 



2 

2.000 
0.625 
0.492 
0.445 

15 
2 

1.999 
0.400 
0.480 
0.519 

13 
2 

1.878 
0.154 
0.260 
0.140 



23 
1 

1.000 
0.000 
0.000 
na 

15 
1 

1.000 
0.000 
0.000 
na 

24 
1 

1.000 
0.000 
0.000 
na 

10 

1 

1.000 
0.000 
0.000 
na 

19 
1 

1.000 
0.000 
0.000 
na 

26 
1 

1.000 
0.000 
0.000 
na 



23 
3 

2.957 
0.739 
0.658 
0.346 

12 

3 

2.943 
0.917 
0.622 
0.127 

19 
4 

3.119 
0.579 
0.633 
0.821 

10 
3 

2.889 
0.600 
0.585 
0.951 

18 
4 

3.594 
0.611 
0.674 
0.338 

12 
3 

2.859 
0.417 
0.538 
0.052 



23 
1 

1.000 
0.000 
0.000 
na 

14 
1 

1.000 
0.000 
0.000 
na 

30 
1 

1.000 
0.000 
0.000 
na 

10 
1 

1.000 
0.000 
0.000 
na 

19 
1 

1.000 
0.000 
0.000 
na 

23 
1 

1.000 
0.000 
0.000 
na 



23 
5 

3.141 
0.478 
0.546 
0.028 

15 
5 

4.239 
1.000 
0.731 
0.017 

30 
7 

3.890 
0.800 
0.687 
0.832 

10 
4 

3.510 
0.500 
0.615 
0.320 

19 
6 

3.572 
0.579 
0.615 
0.546 

23 
6 

3.627 
0.652 
0.601 
0.459 



23 
9 

5.149 
0.826 
0.781 
0.002 

11 

7 

5.612 
1.000 
0.826 
0.001 

11 

6 

5.108 
0.818 
0.793 
0.569 

10 

5 

4.236 
0.700 
0.745 
0.929 

13 
7 

5.174 
0.615 
0.799 
0.279 

17 
4 

3.280 
0.412 
0.538 
0.163 



N is sample size, #A is number of alleles, A R is allelic richness, H D is observed heterozygosity, H E is expected heterozy- 
gosity, P H -w is probability of conformance to Hardy-Weinberg expectation. Sampling localities are shown in Table 1 
and Fig. 1. P-values < 0.05 are in bold. 
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Appendix 2. Pairwise F ST (lower triangle) and P-values (upper triangle) between eight sampling localities of landlocked Atlantic salmon 
and one sample of anadromous Atlantic salmon from the River Namsen watercourse. Significant after corrections for multiple testing. 



Pop 


ANA78 


A 


B 


C 


D 


E 


F 


G 


E55 


E78 


G56 


G/8 


A79 


ANA78 




-0 


~0 


-0 


~0 


~0 


~0 


~0 


~0 


-0 


~0 


-0 


-0 


A 


0.216 




0.61 5 


0.889 


0.891 


0.020 


~0 


~0 


0.001 


0.016 


-0 


-0 


0.940 


B 


0.233 


-0.006 




0.338 


0.090 


0.083 


~0 


-0 


0.000 


0.326 


-0 


-0 


0.227 


C 


0.260 


-0.012 


0.003 




0.713 


~0 


-0 


-0 


0.000 


0.000 


-0 


-0 


0.462 


D 


0.217 


-0.008 


0.026 


-0.003 




0.001 


-0 


-0 


0.001 


0.003 


0.000 


-0 


0.480 


E 


0.228 


0.018 


0.020 


0.040 


0.046 




0.087 


0.021 


-0 


0.360 


0.038 


0.033 


-0 


F 


0.236 


0.060 


0.055 


0.086 


0.082 


0.017 




0.483 


0.018 


0.204 


0.114 


0.052 


~0 


G 


0.254 


0.077 


0.061 


0.095 


0.098 


0.022 


-0.003 




-0 


0.244 


0.106 


0.171 


~0 


E55 


0.225 


0.054 


0.062 


0.073 


0.075 


0.048 


0.019 


0.045 




0.002 


0.001 


~0 


0.000 


E/8 


0.243 


0.035 


0.005 


0.045 


0.060 


0.011 


0.013 


0.011 


0.040 




0.035 


0.076 


0.000 


G 5 6 


0.227 


0.092 


0.078 


0.109 


0.094 


0.036 


0.020 


0.032 


0.051 


0.034 




0.059 


-0 


G78 


0.224 


0.065 


0.049 


0.083 


0.080 


0.019 


0.010 


0.014 


0.045 


0.019 


-0.001 




~0 


A79 


0.251 


-0.010 


0.030 


0.007 


0.008 


0.062 


0.116 


0.137 


0.118 


0.084 


0.151 


0.112 





Sampling localities are shown in Table 1 and Fig. 1. 
P-values < 0.05 are in bold. 
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